r 



1176 00159 9480 


NASA CONTRACTOR REPORT 166302 


A/4S/7 C/e~/66, 


NASA-CR-166302 

19820012373 


Carbon Fibers Conductivity Studies 


Cary Y. Yang and 
Aldona M. Butkus 



P 

r.« 


/l Fi M 

i 


I'M 1 5 1982 

LANGLEY RESEARnv h-.mtfr 
LIBRARY, NASA 

Hampton, Virginia 


CONTRACT NAS2— 1018 8 
August 1980 





NF02348 





NASA CONTRACTOR REPORT 166302 


Carbon Fibers Conductivity Studies 


Cary Y. Yang 

Aldona M. Butkus 

Surface Analytic Research, Inc. 

1662 Parkhills Avenue 

Los Altos, California 94022 


Prepared for 
Ames Research Center 
« ; under Contract NAS2-10188 


NASA 

National Aeronautics and 
Soaro Administration 

Ames Research Center 

Moltett hold. California 94035 







This Page Intentionally Left Blank 



TABLE OF CONTENTS 


Page 

INTRODUCTION ! 

APPENDIX I - Electronic Structure of Clustering Modeling Pan Carbon Fibers . . 6 

APPENDIX II - Properties and Modeling of PAN- Based Carbon Fibers 25 

APPENDIX III - Cluster Model Studies of Oxygen Intercalated Graphite 50 


iii 



INTRODUCTION 


The high tensile strength and light weight characteristics of 
carbon iioer composites have attracted industrial attention as 
potential replacement ior heavier metal materials used in space and 
ground vehicles. However, due to their high electrical conductivity, 
inadvertent release of the carbon fibers can pose fire hazards by 
interfering with unprotected electrical contacts. In order to attempt 
to alleviate potential danger from carbon fibers and yet to preserve 
their properties as light weight, high strength material, a basic 
understanding of the process of electrical conduction in carbon fibers 
and how it is affected by structural variations became desirable. Two 
approaches were used in this study. One focused on the role residual 
nitrogen played in fiber conductivity while the other addressed the 
feasibility of fiber modification to achieve electrical conductivity 
reduction. The first approach supplemented the work by Eabii, e£ 
on electronic energy band calculation for graphite with substitutional 
hetero atoms. The second approach complemented the band structure 
calculations of Fong and Tsang^ on oxygen intercalated graphite. 

When polyacrylonitrile (PAN) is subjected to oxidation and 
v,ub sequent heat treatment, most of the non-carbon elements are driven 
off and the polymer chains condense forming aromatic hexagonal networks 
of carbon atoms. The resultant carbon fiber possesses enormous 
stiffness and strength as well as electrical conductivity several 
orders of magnitude greater than that of the starting material. Recent- 
elemental analysis of PAN samples carbonized between o50-900*C^^ 
indicate substantial nitrogen retention. 3ased on l) the nitrogen 
content results, 2) electrical conductivity as a function of heat- 
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treatment temperature (HTT)^ 3) fiber mean crystallite dimensions ^>6) 
and 4) an idealized scheme previously proposed^' 5 ' for the carboni- 
zation of PM fibers (Figure 1 of Appendix I), nitrogen site models of 
PM fibers at low HIT were studied. 

Calculations were carried out on cluster models of PM fiber 
consisting of carbon, nitrogen, and hydrogen atoms using the modified 
intermediate neglect of differential overlap (METDO) molecular orbital 
(MO) method The models were developed based on the assumption that 
PM fibers with HIT below 1000 ’C retain nitrogen in a graphite-like 
lattice. The clusters studied are shown in Figure 2 of Appendix I. 

For clusters modeling an "edge" nitrogen site (Appendix I, Figure 2, 
clusters a,b,c) analysis of the occupied MO's indicated an electron 
distribution similar to that of graphite. A similar analysis for the 
somewhat less stable "interior "nitrogen site (Figure 2, cluster d) 
revealed a partially localized 'if electron distribution around the 
nitrogen atom. A summary of the electron distributions is presented in 
Figure 3 of Appendix I. The calculated binding energies for the 
clusters are reported in Table II of Appendix I. 

Differences in bonding trends and structural stability between ' 
edge and interior nitrogen cluster calculations modeling PM fiber led 
to the following conclusions. If an interior nitrogen model is adopted 
in order to explain the high nitrogen content in low HTT PM fiber, loss 
of nitrogen from interior sites is predicted to take place preferential 
to nitrogen loss from edge sites. Ho significant change in stiffness 
or crystallite size is expected with this interior nitrogen loss. 
However, this process could be at least partly responsible for the 
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observed dramatic conductivity variation^ in the HTT range of 600- 
1000 C. At HTT above 1000 C the conductivity in crease is leveling off 
and the nitrogen to carbon atomic ratio could be more readily expl ain ed 
by the presence of only edge nitrogens.. A more detailed explanation of 
the eifects of residual nitrogen loss in PAH fibers is presented in 
Appendices I and II. 

The second part of the carbon fiber study dealing with possible 
fioer modification considered the possible electrical conductivity 
reduction oy intercalation as modeled by oxygen intercalated graphite, 
ihis idealized model was chosen based on the observation of extensive 
graphitic domains present in high HTT PAN fibers. Oxygen was chosen as 
the intercalant since it could form strong bonds with carbon and thus 
up the conduction electrons. Electrical conductivity reduction 
would then be expected. 

For i.he oxygen intercalation study, cluster calculations were again 
carried out using the MINDO approach. The graphite models consisted of 
both a full three dimensional representation as well as an independent 
layer model of the lattice. Four different positions were chosen for 
the oxygen atom and are indicated in Figures 1,2, 3, 4 of Appendix III. 

The oxygen-graphite layer distances in the cluster models were optimized 
within a rigid lattice approximation. This results in the oxygen atom 
being the most energetically favored at a distance of 1.25 A above a 
carbon-caroon bond. Additionally , oxygen interaction with the carbons 
ouly one graphite plane predominated. A more detailed account of the 
calculations as well as a summary of the results can be found in 
Appendix HI. Evidence of strong oxygen-carbon bending suggested 
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strong interaction with the graphite conduction electrons, hence 
leading to an insulating phase. 

The calculations that have been performed dealing with the carbon 
fiber electrical conductivity problem represent the first attempt in 
interpreting, elucidating, and predicting properties in carbon fibers 
using a loc ali zed cluster model approach. In view of the general lack 
of long-range order in PM carbon fibers a less extended picture is 
necessary than that provided by band structure calculations and cluster 
models can help establish general bonding trends and local electronic 
charge transfers and binding strengths. Coordination of such theore- 
tical approaches with experimental efforts are essential in attaining 
a better understanding of PM fibers and hence a better control over 
their properties. 
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Appendix I 


ELECTRONIC STRUCTURE OF CLUSTERS MODELING PM CARBON FIBERS 

Aldona M. Butkus and Cary Y. Yang 

Surface Analytic Research, Inc. 

1662 Parkhills Avenue 
Los Altos, CA 9 U 022 

ABSTRACT 

In an attempt to understand the process of electrical conduction in 
polyacrylonitrile (PAN) -based carbon fibers, calculations were carried 
out on cluster models of the fiber consisting of carbon, nitrogen and 
hydrogen atoms using the modified intermediate neglect of differential 
overlap (M3ND0) molecular orbital (MO) method. The models were 
developed based on the assumption that PM carbon fibers obtained with 
heat treatment temperatures (HIT) below 1000 ’C retain nitrogen in a 
graphite-like lattice. For clusters modeling an "edge" nitrogen site, 
analysis of the occupied MO's indicated an electron distribution similar 
to that of graphite. A similar analysis for the somewhat less stable 
"interior" nitrogen site revealed a partially localized IT electron 
distribution around the nitrogen atom. The differences in bonding 
trends and structural stability. between edge and interior nitrogen 
clusters led to a two-step process proposed for nitrogen evolution with 
increasing HTT. 
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I. INTRODUCTION 


When polyacrylonitrile (RAN), (CH 2 -CH(CN) ) n where n is typically 2000 
in commercial fibers, is subjected to oxidation and subsequent heat 
treatment, most of the non-carbon elements are driven off and the 
polymer chains condense forming aromatic hexagonal networks of carbon 
atoms. The resultant carbon fiber possesses enormous stiffness and 
strength as well as electrical conductivity several orders of magnitude 
greater than that of the starting material. Electron microscopy and 
x-ray diffraction studies (1) - of these carbon fibers indicate that each 
fiber consists of ribbons of crystallites bundled together and cross- 
linked, while being aligned along the fiber axis. The typical crystal 
structure within each ribbon is that of turbostratic graphite 
possessing no regular stacking order in the direction perpendicular to 
the basal plane. Since PAN carbon fibers are non-graphitizable, the 
interplanar spacing, while decreasing with increasing heat treatment 
temperature (HIT), never reaches that of pure graphite even at the 
highest HTTs. Values of and L c , the mean crystallite dimensions 
parallel and perpendicular to the basal planes, increase with Increasing 
HTT. For PAN carbon fibers with HIT around 1000*C, values of about 20A 
for the width and about 1QA for the thickness are typical. 


Recent elemental analysis of PAN samples carbonized between 650-900* C 
indicated substantial nitrogen retention^) the nitrogen to carbon 
atomic ratio ranging from 0.24 to 0.15. An idealized scheme showing the 

nitrogen sites previously proposed (4 ' 5) far the carbonization process of 
PAN fiber is indicated in Figure 1. Assuming the ribbon models shown in 
Figure 1 as possible intermediate structures during the 
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carbonization process, we have calculated the nitrogen to carbon ratio 
with increasing ribbon width, L ft , and the results are given in Table I. 
This model allows the presence of only "edge" nitrogens. In order to 
accomodate the high percentage of nitrogen in the sanples reported in 
Reference 3, L a have to be between 10 and ISA in that range of HIT. 
However, the crystallite dimensions of these samples are not available. 
Since extensive graphitized domains exist in PAN carbon fibers with HIT 
below 1100 C , such small widths seem unlikely. Hence the retained 
nitrogen probably resides in additional sites. Discounting the presence 
of interstitial nitrogens, one is left with the only alternative, namely, 
the substitution of an interior carbon in the two-dimensional graphite 
environment by a nitrogen atom. This substitution is similar to that in 
heterocyclic compounds, with bridgehead nitrogen atais, e.g. 
cycl [3.3.3] a z ine s . However, the only, chemical evidence for the 
existence of a stable interior nitrogen during the carbonization pr ocess is 
the large amount of nitrogen at lew htt. The effects on structure and 
electrical conductivity of the presence of such, an interior nitrogen in 
contrast to an edge nitrogen constitute the main thrust of this work. 

The role of other mechanisms such as hydrogen loss or contribution from 
impurity levels ' giving rise to changes in electrical conductivity 
during the carbonization process of carbon fibers has not been address ed 
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II. METHOD OF CALCULATION 


In order to study the structural stability and general bonding trend of 

nitrogen substitution in PAN carbon fibers, modified intermediate 

neglect of differential overlap (MINDO/3)^ molecular orbital 

calculations have been carried out on clusters of atoms consisting of 

carbon, nitrogen and hydrogen modeling the nitrogen site environment. 

The MBIDO/3 approach was chosen because it allows one to obtain the 

ground state solution tq the electronic Schrodinger equation for 

polyatomic systems in considerably less computation time. This savings 

over the ab initio Hartree-Fock-Roothaan molecular orbital approach is 

achieved through either the neglect or the replacement by parameters of 

certain Hartree-Fock Hamiltonian matrix elements. These parameters have 

been previously determined through a lengthy iterative procedure by 

fitting experimental geometries and heats of formation for a large class 

( 8 ) 

of molecules . Additionally, since it is of interest in the present 
study to deteimine possible changes in geometry with the introduction 
of nitrogen into the graphite lattice, total energy minimization as a 
function of atomic coordinates would be desirable and could be economically 
performed with the MINDO/3 approximation. 

The clusters chosen to model the nitrogen sites in the fiber are shown 
In Figure 2. A few calculations were carried out on three dimensional 
cluster models of the nitrogen site assuming an ideal graphite stacking 
geometry. But since these three dimensional studies yielded results 
identical to those of the two dimensional models and the interlayer 


9 



structure of the fiber is basically turbootratic , with no long range 
order in the direction perpendicular to the basal plane , the ensuing 
analysis will be based on the two dimensional models. In order to 
minimize edge effects introduced by a finite cluster representation of 
the lattice, hydrogen atoms were used to saturate the valence of the 
outer carbons. Geometry optimizations were carried out and the heats 
of formation were obtained. Study of the effect of nitrogen 
substitution on conductivity was based on the calculated electron 
distribution in orbital3 perpendicular to the basal plane (IT electrons),. 
This was accctnplis^bd using a Mulliken population analysis^) which is a ' 
means of partitioning electron density among the various atoms in the 
cluster from the calculated molecular orbitals. 

Previously, theoretical study of the one-ring and three-ring structures 
had been carried out. The single ring calculation^^ was 'performed with 

MINDO/3, which is also employed here. But since the Mulliken 

( 10 ) 

population analysis was not reported , we have repeated the one-ring 
calculation. The earlier study of the three-ring structures was done 
using a more approximate semi-empirical approach with only the 7T 
orbitals determined self-consistently^ ^ • 
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III. RESULTS AND DISCUSSION 


Geometry optimization of the various cluster models shown in Figure 2 
yields planar environments about the nitrogen for both the interior and 
the edge nitrogen models. The carbon -nitrogen bond length for the edge 
nitrogen model varies from 1.335A (single-ring model) to 1 . 356 a (four- 
ring model), while a carbon-nitrogen bond length of 1.4l6A is obtained 
for the interior nitrogen case. The carbon -hydrogen distances are all 
approximately 1.1QA, a value typically found in aromatic molecules. 
Differences in geometry between the various edge nitrogen models are 
within 2$ indicating that any structural changes due to the presence of 
nitrogen are local to the substitution site and do not significantly 
affect the adjacent rings. 

Analysis of the occupied molecular orbitals for the edge nitrogen models 
shows that as in the pure graphite case one electron from each ring atom 
(carbon or nitrogen) is present in the it -type orbitals perpendicular 
to the basal plane, with the remaining electrons participating in 
O' -bonding with the adjacent in-plane carbons. The non-bonding lone 
pair orbital on nitrogen normal to the edge surface then becomes the 
analog of the third O' bond of the pure graphite case. 

For the interior nitrogen site, three O' - type bonds are formed with 
the adjacent carbons in the basal plane similar to those in graphite, 
the carbon-nitrogen bond distance being very close to that of carbon - 
qarbon in graphite. However, the remaining two valence electrons, 
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being partially localized on nitrogen, are donated to the IT system, 
resulting in a very different 4 Tt electron distribution from that of 
graphite. In order to better see these changes in the It system, the 
fractional If electron charge at each atom is calculated by dividing 
the electron population in the orbital perpendicular to the basal plane 
(p ff ) by the total atomic population for the particular atom. A 
summary of the calculated values is given in Figure 3* This fractional 
It charge at each atom is used instead of the explicit p ff population 
because the former approach better represents the actual d'-IT 
breakdown for all atoms, regardless of their environment. 

In order to investigate the relative stability of the nitrogen in the 
two sites, the average binding energy per atom relative to atomization 
is calculated for an edge nitrogen model (4-ring) , an interior model 
(3- ring), as well as a pure graphite model analogous to the edge model. 
The results are summarized in Table II for the three cases. 

From the fractional IT electron distribution shown in Figure 3 , it 
becomes apparent that the edge nitrogen gives rise to a more even 
electron distribution than the interior nitrogen case. In fact, the 
former resembles pure graphite in its TT electron distribution. As 
with the geometry optimization results, effects due to the presence of 
an edge nitrogen are confined to the nearest carbon neighbors. In 
contrast, the Tt electron distribution for the interior nitrogen site 
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shown a significant variation an one passes from atom to atom. iJinco 
it is the electrons in the TT system that are responsible for the 
conductivity along the basal plane, the partial localization of the 
electrons would imply a reduction in conductivity for the fiber with an 
interior nitrogen. On the other hand, the edge nitrogen would not be 
expected to play as large a role in changing the conductivity. 

The relative stability of the two sites as estimated from the average 
binding energies yields an energy difference of approximately 0.25eV 
(5*8 kcals/mole) between the edge and interior models indicating that 
the edge site is energetically favored over the interior one. This 
energy difference is remarkably close to the experimentally determined 
value of 5*9 kcals/mole for the barrier to inversion for ammonia. The 
transition state for this process requires the atems to be coplunar, 
placing a steric constraint on the nitrogen similar to that present in 
the interior nitrogen model. 

Thus if an interior nitrogen model i3 adopted in order to explain the 

high nitrogen content in low HTT PAN fiber, loss of nitrogen from 

interior sites is predicted to take place preferential to nitrogen loss 

from the edge sites. No significant change in stiffness or- crystallite 

size is expected with this interior nitrogen loss. However, we propose 

that this process is at least partly responsible for the observed 

(\o) 

dramatic conductivity variation with heat treatment of TAN fibers 
iu the temperature range of 600-1000*0. If one adopts the two- 
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dijnenoional graphite touching band picture ^* 3 \ the partial localization 
of 7T electrons most likely results in a somewhat flattening of the 
valence 'and conduction bands. According to semiconductor band theory, 
the consequences of such band flattening are generally two-fold. 

First, the band gap opens up; and second, the carrier effective masses 
increase. One would also expect that nitrogen in the graphite lattice, 
whether or not in an ordered arrangement, will shorten the carrier 
life-time» through increased charge localization*, Hence, the overall 
effect of nitrogen retention on the electrical properties of FAN carbon 
fibers is a decrease in the conductivity resulting from simultaneous 
decreases in carrier densities and mobilities. In view of the more 
localized TT electrons around the interior nitrogen atom, we believe 
that nitrogen retention in the interior provides one important physical 
mechanism responsible for the drastic decrease in conductivity of FAN 
carbon fibers up to 1000 *C HIT. 

. (u) 

At HIT above 1000 C the conductivity increase is leveling off and 

. t m ) 

values for L (20-3QA) and the nitrogen to carbon atomic ratio' 

3 » 

(0.09) can now be readily accounted for by the presence of only edge 

v 

nitrogens. The results reported here do not indicate a direct effect on 
the conductivity due to the loss of edge nitrogen. However, it can be 
inferred from the largo concentration of electron density about the 
edge nitrogen atom (due to the protruding lone pair) that coalescence 

1 • 

of the ribbons to form a larger hexagonal carbon net would be inhibited 
based on electrostatic grounds. With the loss of these edge nitrogens 
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leading to coalescence and hence larger crystallite dim ensions, the 

carrier mobility due to ribbon boundary scattering will increase. 

Since in this range of HIT, the resistivity is probably dominated by 

( 15) 

scattering from ribbon boundaries , and since larger crystallite 
dimensions give rise to greater ft delocalization as well as extended 
O' linkage, increases in conductivity and stiffness would result. 

This proposed two-step nitrogen elimination could also help explain the 
observed sharp increase in nitrogen evolution^ around 900 *C, since 
at this temperature nitrogen loss from both sites could be occurring 
simultaneously. 
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I V. CONCLUSIONS 


The calculations presented here represent the first attempt in 
interpreting and elucidating properties in carbon fibers using a 

(n is) 

localized cluster model approach. Previous attempts * have been 
based on information directly derived from an assumed band structure. 

It is our contention that a less extended picture is necessary to 
provide complementary data whereby properties of the material, can be 
better understood. In view of the general lack of long-range order in 
PAN carbon fibers, cluster models can help establish general bonding 
trends such as those found in the nitrogen-containing clusters. 

Our suggestion that the If bands are flattened due to nitrogen 
retention will hopefully be verified by the ongoing calculations of 

/ IQ \ 

Rabii on nitrogen graphite. From the band structure, we hope to 
est imat e the conductivity as a function of nitrogen concentration, i.e. 
HTT. The two-step nitrogen elimination model proposed here should 
stimulate more measurements of the crystallite dimensions of low HTT 
PAN carbon fibers, as well as more spectroscopy experiments. We 
believe that systematic experimental and theoretical efforts, properly 
coordinated, are essential in attaining a better understanding of PAN 
carbon fibers, and hence a better control over their properties. 
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TABLE 1. Nitrogen to / c;irbon atomic ratio and number of interior ring:: 
as a function of approximate ribbon width L for ribbon model 
of PAN fiber carbonization process. 


\0i) 



Ratio N to C 

Number interior rln^c 

7 

0.33 

i 

9 

0.25 

2 

11 

0.20 

3 

13 

0.17 

h 

15 

0.14 

5 • 

20 

o.u 

7 

30 

0.07 

12 
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TABLE II. Average binding energy per atom for the nitrogen edge 

(4-ring) and interior (3-ring) models and the four- ring 
cluster graphite model. 


Cluster 
N-edge 
N- interior 
Graphite 


Energy (eV/ atom) 
5.28 
5.03 
5.22 
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FIGURE CAPTIONS 


Figure 1. 
Figure 2. 
Figure 3* 


The development of graphite structures during the pyrolysis 
of thermally stabilized PM. 

Cluster models chosen for edge (a,b,c) and interior (d) 
nitrogen sites in PAN fiber. 

Fractional ir electron polulation for (a) edge nitrogen 
model, (b) interior nitrogen model, and (c) graphite model. 
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(d) 

Figure 2 
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Figure 3 
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Summary 

A review of the properties of PAN-based carbon fibers is presented 
with a view toward developing a meaningful model for theoretical calcu- 
lations and the analysis of experimental data. After reviewing the 
status of past and present models, the results of the first cluster Cal- 
culations of the electronic structure of carbon fibers is presented. 
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Introduction 


There has recently been increased interest in the electrical properties 
of carbon fibers. Much of this interest stems from the highly desirable mechani 
cal properties of these fibers (four times stronger than steel and only one- 
quarter the density) which leads to their use as strength-giving elements in 
plastic composite materials. Current applications of carbon— fiber composites 
range from tennis rackets to the bay doors on NASA's space shuttle. The possi- 
bility of making strong light-weight cables and other electrical, applications 
provides strong impetus to the materials engineer and solid-state scientist 
to understand the nature of the electronic conduction mechanism in carbon 
fibers and the relationships between the conductivity, the structure, and the 
mechanical properties of the fibers^ all of which are related to properties of 
the precursor material and the processing history. 

The room-temperature electrical conductivities of commercial carbon fibers 
range from 500 - 5000 (Ocm)^ - . These values are considerably lower than the 
basal-plane conductivity of pristing graphite or copper, which have values of 
25,000 and 580,000 (Q-cro)’^, respectively; however, certain chemical treatments 

J 

of carbon fibers have been shown to increase the fiber conductivity by an order 
of magnitude [1-4], into a range with more engineering potential. The fact 
;that this increase- can-' only be produced with fibers which have been processed 
in a certain manner gives further reason to study the basic electrical proper- 
ties of carbon fibers. 

In this paper we first present a review of the properties of carbon fibers 
with a view toward developing a meaningful model for theoretical calculations 
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and the analysis of experimental data, we then discuss the status of past 

and present models and present the results of cluster calculations which provide 

the first theoretical picture of the electronic structure of carbon fibers. 

For the sake of brevity and in order to develop a specific model, we have 
chosen to concentrate on carbon fibers made from precursor fibers of poly- 
acrylonitrile (PAN). PAN-based carbon fibers are presently used for most com- 
posite applications, having replaced the rayon-based carbon fibers which were 
developed a decade ago. Pitch-based fibers, which are still in the develop- 
mental stage, are the third type of carbon fiber of significant commercial 
interest. Even though we concentrate on the PAN— based fiber many aspects of 
the discussion will apply to carbon fibers from other precursors, as will be 
pointed out. 

2. Properties of PAN-Based Carbon Fibers 

2.1 Production, Chemistry, and Mechanical Properties 

The homopolymer polyacrylonitrile (PANJ upon which the carbon fiber' is 
largely based, has the formula [CH2 _ CH(CN) ] , where the degree of polymeri— 

* 

zation, n, for commercial fibers is usually higher than 2000. Commercial PAN 

t 

fibers such as DuPont f s "Orion" also contain copolymers such as vinyl acetate, 
vinyl chloride, or styrene to improve such properties as their ability to 
stretch or their affinities for dyes, and the precursor materials for commercial 
carbon fibers also contain a number of additives [5]. 

Commercial carbon fibers are produced from PAN fibers by first "oxidizing" 
the PAN fiber under tension at a temperature near 250°C, then carbonizing and 
heat treating in an inert atmosphere at temperatures from 1200 to 3000°C. 
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Figure 1 shows some possible intermediate structures during these production 
steps as conceived by Watt [6]. The lower temperature "oxidation" drives off 
hydrogen, possibly produced ketonic groups (-C0-) , and causes a cross-linking 
of CN groups to produce a ladder structure. Tension is applied to assure 
alignment of these ladders along the fiber axis, a necessary condition for the 
production of high modulus fibers. This alignment is preserved with the final 
heat treatment even if the tension is removed. Between the 250° C "oxidation" 
and the subsequent high temperatures treatments various amounts -of water, HCN, 
N 2 , CO, CO 2 and other gases are given off and the ladder chains evolve into an 
aromatic hexagonal network of carbon atoms. Figure 2, from Fiedler, et al. 

[7] shows typical data for the rate of gas evolution for PAN fibers as various 
heat treatment temperatures (HTT) . Fibers produced with HTTs near 1500°C 
typically have from 5 to 16 percent non-carbon residue. The exact structure 
of the carbon sheets is not known and the way that the residue fits into 
the actual atomic structure of the fiber has not yet been clearly determined. 
One of the goals of our research is to gain an understanding of this structure 
through a correlation of experimental results with theoretical models. 

j 

The processing HTT can be directly correlated with the resulting mechani- 
cal properties of the fibers as shown by the data of Watt and Johnson [8] , pre- 
sented in Figure 3. The modulus of PAN-based carbon fibers rises continuously 
with HTT, however, the tensile strength passes through a maximum value near 
1500°C and then decreases with increasing HTT. Because of this maximum, 
commercial PAN-based carbon fibers are available either with high tensile 
strengths and medium stiffness moduli, or high moduli and medium strengths. 

The former group includes such fibers as Union Carbide’s Thomel 300, manufac- 
tured by Toray in Japan; Celanese’s Celion 6000, from Japan by Toho-Beslon; 
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Figure 1. 
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The chemical changes which occur as PAN fibers are turned into 
carbon fiber as described by Watt, Ref. 6. 
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Figure 3 The HTT dependence of tensile strength and Young's modulus for 
PAN-based carbon fibers, from Ref. 8. 
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and Hercules' HTS2 until recently produced in England by Courtaulds. An 
example of a high modulus fiber is Celanese's GY70. Table 1 summarizes 
typical values for moduli and tensile strengths and their relationships with 
the HTT along with several other properties of PAN-based carbon fibers and the 
properties of graphite. Whereas PAN-based carbon fibers cannot be trans- 
formed into graphite it is interesting to see how the density, modulus, and 
resistivity of the fibers approach those of graphite as the HTT is increased. 

This relationship will be discussed further in subsequent sections. 

2.2 Structure 

Typical of the models of the structure of PAN-based fibers is that pre- 
sented by Wicks [9], shown in Figure 4, in which a fiber consists of bundles 
of microfibrils that lie along the fiber axis. In this model, the result of 
numerous x-ray and electron diffraction studies, each fibril is made up of 
parallel layers of hexagonally-arrayed carbon atoms. The layers are "turbostratic" 
i.e. randomly stacked ; and lack the three-dimensional order of graphite as 
indicated by the absence of certain (hkl) x-ray reflections. The details of 
the structure of each fibril, as well as the alignment of the fibrils, depend 
upon the HTT and the stress conditions during production. For example, the , 
separation between the carbon sheets varies from 3.4 to 3.7 angstroms [10], with 
the closest spacing for PAN fibers reached only with heat treatment above 
2000°C. However, even at the highest treatment temperatures the ideal graphite 
spacing, 3.35 angstroms, cannot be attained. This is another way of saying that 
PAN fibers are non-graphitizable as we mentioned earlier. Fibril thicknesses, 

L , vary from 30 to 100 angstroms, with a characteristic straight section, L , 
of the order of 100 angstroms, again, depending on the heat treatment history. 

A typical set of data for L and L from Johnson and Tyson [11] is plotted in 

c a 

Figure 5. 
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In Figure 4, needle-shaped voids can be seen. These voids, which are 
mostly inaccessible to helium, can represent from 5 to 30 percent of the total 
volume. The fibers also contain impurities and residuals from the precursor 
materials, e.g., the high strength, medium-modulus PAN-based carbon fibers 
are only 90 percent carbon, as mentioned earlier. 

Recently, several heterogeneous structures have been proposed for PAN- 
based carbon fibers, such as a core with a radial structure of crystalline 
webs separated by large voids and the core surrounded by a highly crystalline 
sheath [12]. However, Bennett and Johnson [13] have shown that these fibers 
have a rather homogeneous structure not unlike that depicted in Figure 4. 

Only with heat treatment at 2500°C is the sheath/core heterogeniety observed 
and then the sheath is thin, from 100 to 1500 angstroms. 

The modulus and strength of carbon fibers is basically derived from the 
covalent bonds between the hexagonally-arrayed carbon atoms within each sheet. 
These are the same bonds which give graphite a basal-plane modules or 150 X 10 
psi. The planar covalent bonds, called a - bonds, are sp hybrid orbitals 
involving three of the four valence electrons of each carbon atom. They are 
actually stronger than the sp bonding orbitals of diamond. In order for these 
- bonds to be effective the sheets of carbon atoms must be parallel to the fiber 

i 

axis. This is achieved by stretching the fiber during the initial oxidation 
assuring that theiladder chains are aligned along the fiber axis. 

2.3 Electrical Properties 

The relationship between the HTT and tjie electrical properties of PAN-based 
carbon fibers was first published by Yamaguchi in 1963 [14]. His results and 
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a comparison with soft carbons are shown in Figure 6. The dramatic decrease 
in the resistivity above 600°C parallels the onset of the planar— sp^ (or ct) 
bonding which eventually leads to the turbostratic structure which we have 
previously discussed. This relationship between bonding and electrical con- 
duction was clearly shown by Ezekial in 1970 [15] when he comapred the resis- 
tivity and Young's modulus for both rayon-based and PAN-based carbon fibers. 
Ezekial s results are shown in Figure 7 which also includes data points from 
Table 1. The information from Figures 6 and 7 leads us to picture the increased 
conduction as arising from the formation of tt — bands which accompany the 
formation of the strength-giving a bonds. Pristine graphite, the ideal crys- 
talline material toward which the fibers tend, is a semimetal with electron and 
hole concentrations of n = p ~ 10^ cm ^ , each with mobilities near 10^ cm^/V-s, 
resulting in a room-temperature resistivity of p = 4 X 10“^ £J-cm. The 

minimum resistivity for the fibers shown in Figure 6 is approximately twenty 
times that value. 


At two recent meetings the first comprehensive data was presented on the 
temperature dependence of electrical conductivity of PAN-based carbon fibers. 
Kalnin and Goldberg [4] gave data for five different sets of PAN-based fiber's, 
which had been .teat treated above 2000°C and had moduli between 50 X 10^ 
and 110 X 10 psi. Gillespie, Ehrlich, and Waters [3] presented p(T) curves 
for the commercial fibers AS, HMS, T-300, and GY-70 as well as two pitch 
fibers. Between 50 and 300 K all of the resistivity curves show a monotonic 
decrease in resistivity with increasing temperature similar in slope and mag- 
nitude to the data resistivity taken by Bright and Singer [16] for pitch-based 
carbon fibers, which is given in Figure 8. The data for the AS and T-300 
fibers, which have HTTs near 1500°C, have the extra feature that p peaks near 
25K with a 4 percent decrease in p between 25K and 4K. The decrease in 
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Figure 6 


Carbon resistivity as a function of HTT, from Ref. 14 
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Figure 7 / 
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The relationship between resistivity and modulus for two types 

* 

of carbon fibers, from Ref. 15. 
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Figure 8 Resistivity as a function of temperature for pitch-based carbon 
* fibers, from Ref. 16. 



resistivity with increasing T most likely indicates that the charge carriers 
in the fibers are thermally activated as in a semiconductor or semimetal with 
shallow overlap of bands. A similar temperature dependence was noted by Hdrinckx 
for rayon-based carbon fibers [17]. The only other published resistivity data 
for PAN-based carbon fibers is that of Robson et al. [18], shown in Figure 9, 
where the resistivity ratio between liquid nitrogen and room temperature is 
given as a function of HTT. There is an interesting similarity between the 
data in Figures 8 and 9: they both have a minimum in the resistivity ratio 
near HTT 1750°C. 

The magnetoresistance effect has proven to be a useful analytical tool in 
semiconductors because its magnitude is related to the mobilities of the carriers. 
Unfortunately, the interpretation of the magnetoresistance in carbon fibers 
and pyrocarbons can be complicated by a negative contribution which is still 
not understood even though it has been studied extensively [19, 20]. Robson 
et al. [18] showd that the negative ef fect&rises in PAN/carbon fibers which 
have HTTs greater than 1750°C, as shown in Figure 10. This is a property 
that will have to be addressed in future models of these materials. 

* 

3. Models , 

Akamatu and Inokuchi [21] showed that compounds with carbon-ring structures 
exhibited conductivities with seraiconductor-like temperature dependences even 
for as few as nine rings, and Mrozowski [22] extended this idea to explain the 
experimental results on soft carbons. He proposed a variable-band-gap model 
in which successively higher heat treatment of soft carbon change them from 
semiconductors to semimetals, with the band gap decreasing to zero and then 
into an overlap region. In Mrozowski T s model the carbon atoms begin to form 
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Figure 10 


The magnetoresistance for PAN-based Carbon fibers as a function 
of HIT, from Ref. 18. 
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an aromatic structure upon heat . treatment near 600°C, with the rr electrons 
forming valence and conduction band with an energy gap of ~0.3 eV. With 
higher heat treatment the gap decreases, going to zero near 1200°C, as the 
lamellar carbon planes become larger and the spacing between the turbostratically 
aligned layers decreases* With heat treatment above 3000°C the carbons assume 
the three-dimensional ordering of graphite, where the bands overlay by approxi- 
mately 0.04 eV. 

Over the years, a number of models have been proposed for carbon fibers 
in attempts to explain and elucidate various structural and electronic pro- 
perties. Yazawa [23], for example, employed a two-dimensional band model of 
graphite with adjustable parameters and successfully predicted the variation of 
raagnetoresistance and Hall coefficient with HTT, magnetic field and sample 
orientation. Both models are based on the energy band picture of a crystalline 
solid. In view of the general lack of long-range order, even for high HTT 
carbon fibers, a less extended picture fbrvthe material needs to be developed, 
not necessarily as an independent verification of the band models, but rather 
to complement them. It was in this spirit that calculations were initiated 

* 

in 1979 by Butkus and Yang on the electronic structure of jmodel clusters of 
PAN carbon fibers, and by Rabii on the band structure of nitrogen graphite. 

9 

The cluster calculations of Butkus and Yang [24] have been performed 
using a combination of two semi-empirical quantum chemical methods, MINDO/3 
and HAM/3. The object of these calculations is to examine the localized bonding 
picutre in nitrogen-retained PAN carbon fibers. Previously, we have established 
i that regions of turbostratic graphite do exist within PAN-based carbon fibers 
at low HTT*s. This gives rise to the underlying assumption that PAN-based 
carbon fibers consist of crystallites with the exact two-dimensional graphite 


44 



ATOMIC CLUSTER MODELS OF 
CARBON FIBER ELECTRON STRUCTURE 




INTERIOR N SITE 

2 p. ELECTRONS HIGHLY LOCALIZED 
3 (O BONO ELECTRONS 
NON-BOCKLEO RINGS 


EDGE N SITE 

1 p. ELECTRON SLIGHTLY LOCALIZED 
ZipfiONO ELECTRONS 
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Figure 11 Two different 


nitrogen locations used for the cluster calculations 



structure, where nitrogen is retained through direct substitutions in the 
lattice. The resulting model clusters automatically fall into two categories. 

One can model either an "interior” nitrogen, or an "edge" site, as shown in Fig. 11. 
Geometry optimization of the clusters yields planar environments for both sites. 

The edge site, however, has a slightly higher average binding energy. An 
analysis of the electron distributions of the valence orbital reveals that the 
edge nitrogen bonding to • its carbons is similar to that of an edge carbon-hydrogen 
pair. However, the interior nitrogen, while maintaining similar bonding on the 
basal plane, results in rather localized n electron distribution around itself. 

To explain these cluster findings in the context of an extended band pic- 
ture would seem presumptous at best. However, while the band calculations of 
Rabii [25] are still in progress, we regard for the time being the conclusion 
of Butkus and Yang on electrical conductivity of PAN carbon fibers as one 
plausible qualitative explanation. 

Based on these results, and the well-known two-dimensinal graphite band 
structure, Butkus and Yang [24] suggest a flattening of the tt bands due to \ess 
tt delocalization. The consequences, in the somewhat simplified picture of semi- 
conductor band theory * are two fold. First, the effective masses increase, 
decreasing the carrier mobilities. Second, the tt - band flattening inevitably 
gives rise to a larger band gap, resulting in lower carrier densities. Hence, 
nitrogen in carbon fibers, especially the interior one, constitutes one physical 
mechanism by which conductivity is reduced. 
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4. Conclusions and Discussion 

It is clear that we are presently only at the beginning of our under- 
standing of the electronic properties of carbon fibers. We have gained a 
mass of information of the structural properties of these materials, mainly 
because of the commercial interest in their use in composites, but we lack • 
a definitive understanding of their basic electronic structure. 

Basic to the problem is the lack of crystalline order. In one sense the 
electronic response to temperature variation is similar to that of a semi- 
conductor, yet on the other hand, we know that increasing conductivity with 
increasing temperature is also a property of an amorphous material. How 
can we discern the true nature of these materials? 

The attempt here has been to review the basic knowledge of the properties 
of PAN-based carbon fibers and to couch these properties in the language of 
simple models. The cluster calculation is the first fundamental theoretical 
investigation into the nature of PAN-based carbon fibers. We hope that it 

shed light for coming investigators. More research is definitely needed. 
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Appendix HI . .. 

Cluster Model Studies of Oxygen Intercalated Graphite 

Aldona M. Butkus 

Surface Analytic Research, Inc. 

1662 Parkhills Avenue 
Los Altos, CA 94022 

Abstract 

Electronic structure calculations were carried out on clusters of 
atoms modeling oxygen intercalation of graphite using the modified 
intermediate neglect of differential, overlap (MINDO) molecular 
orbital method. The graphite models consisted of both a full three 
dimensional representation as well as an independent layer model of 
the lattice. Pour different positions for the oxygen atom were 
considered and consisted of 1 ) an oxygen directly above a carbon of 
one layer and directly below that of the second (site A), 2) an 
oxygen above a carbon of one layer and below the center of a 
hexagonal ring of the second layer (site B), 3) an oxygen above a 
carbon-carbon bond and below a position midway between second nearest 
neighbors of the hexagonal ring (site 0), and 4) an oxygen above 
and below a point midway between second nearest neighbors of the 
hexagonal rings in both layers (site D). The oxygen-graphite layer 
distances in the cluster models were optimized within a rigid lattice 
approximation resulting in the bond site being the most energetically 
favored position for the oxygen atom at a distance of 1 .25A above 
the carbon-carbon bond. Additionally, oxygen interaction with the 
carbons of only one graphite plane predominated not onJ y for the 
most favored position but for all four sites considered, with the 
three dimensional results differing little from those obtained 
using the independent layer model for the various sites. 
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I . Introduction. 


Graphite consists of planes of linked hexagons of carbon atoms 
parallel to one another along the principal C,axls. Under certain 
conditions it is able to take up atoms, ions, or molecules with 
the lattice structure largely unchanged. These intercalation 
compounds have been the focus of both experimental and theoretical 
interest not only for their unusual electronic and chemical properties 
but also because of their technological possibilities. 1 The work 
reported here- is a theoretical study of the electronic structure 
of the oxygen graphite system. 

The graphite orystal structureis well known, consisting of 
hexagonal unit cells containing four atoms. The stacking sequence 
is of ABAB... type where alternate layers are equivalent. Within 
the layer planes, each carbon atom is surrounded by three other 
carbons at a distance of 1 .42A. The separation between the layers 
is 3.35A and thus clearly indicates that interaction between atoms 
in adjacent planes is less pronounced than between intra-plane 
carbons. If in-plane as well as off -plane nearest neighbor atoms 
are considered, the graphite crystal structure gives rise to two 
different types of atomic environments. In one type of environment 
an atom has three in-plane neighbor atoms and carbons directly 
above and below it in adjacent planes. The other type of atom site 
does not possess the latter two atoms. 

In contrast, little is definitively known about the structure 
of the oxygen Intercalation compounds of graphite. The reported 
carbon/oxygen ratios are never lower than two and range from 2.4 
to 4.0 depending on the degree of oxidation. 2 X-ray studies indicate 
that the separation of carbon layers falls within the range of 
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6 - 11 A. 3 There la some controversy regarding the planarity of the 
carbon layers since in view of the absence of metallic properties 
it is believed that the layers would not be planar as in graphite. 
The carbon atom now attached to four ligands would be tetrahedral 
causing the layers to buckle. However, more recent electron 
microscopy and x-ray diffraction studies done by Carr 4 suggest 
that the layers in oxygen graphite are flat.. 

For the present cluster model studies of oxygen intercalated 
graphite, the graphite layers were assumed to be flat and the atoms 
at their equilibrium positions. Attention was directed toward the 
determination of the most energetically favored site for the oxygen 
atom and the elucidation of the type of bondin'? without any lattice 
relaxation allowed to take place. Oxygen intercalation in the dilute 
limit was modeled since interaction between oxygen atoms was not 
considered. 

II . Method of Calculation 

The present study of the electronic structure of oxygen 
Intercalated graphite consists of using clusters of carton, hydrogen 
and oxygen atoms to model the oxygen-graphite system and carrying 
out modified intermediate neglect of differential overlap (MINDO/3) 
molecular orbital calculations. 5 This particular molecular orbital 
approach, though semi -empirical in nature affords the ability to 
study large molecular models of the graphite lattice with sufficient 
economy, providing a useful approximation to the solution of the 
Ilartree-Fock equations. The savings in computational requirements 
:1s achieved by the replacement of certain difficult integrals by 
appropriately determined parameters. The particular parameterization 
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in the MIN]X)/3 approach reliefs on the accurate reproduction of 
experimental geometries and heat of formation of molecules 
containing the appropriate molecular species.® 

In choosing the cluster models of the oxygen intercalanl. sites, 
preliminary calculations were carried out using the Independent layer 
model of the graphite lattice in order to determine possible 
oxygen binding sites. Models for the graphite layer ranged from a 
single ring of six carbon atoms to a condensed four ring system. 
Hydrogen atoms at a carbon-hydrogen distance of 1 . 1 I were used to 
saturate the dangling o' carbon bonds at the edges of the cluster 
models. Such a treatment of edge carbons was net expected to introduce 
great error since the oxygen-lattice interaction should mainly take 
place with the IT electron system. Additionally, energy differences 
between the lattice model and the lattice-oxygen intercalant model 
were of interest so that cancellation of errors due to the edge 
treatment was expected to occur. Use of the independent layer model 
of the lattice further allowed the determination of the minimum 
cluster size necessary to adequately^ repi’ssent a property of interest. 
Comparison of the present two dimensional lattice MIHDO results 
with the only theoretical studies of the oxygen-graphite system done 
previously concerned with chemisorption^’® ’ 9 was a i so possible. 

Pour different positions for the oxygen atom in the graphite 
lattice were considered and the cluster models chosen for study 
■are indicated in Figures 1,2,3 and 4. The first site, site A, places 
the oxygen directly above and below carbon atoms. The second, site B, 
has the oxygen directly above a carbon and below the center of a 
hexagonal ring of the second layer. Thus in the independent layer 
models of this site, clusters are considered with t.ho oxygen above a 

\ 


53 



5 


carbon aa well ns .above the center of a ring of carbons separately. 

The oxygen can also be positioned above a carbon-carbon bond, site C, 
which results in a position below a point midway between second 
nearest neighbors of the second ring in the three dimensional lattice 
model. The last oxygen-graphite site considered consists of the 
oxygen above and below the point midway between second nearest 
neighbors of the hexagonal rings. 

For all of the oxygen-graphite sites considered, MINDO calculations 
were first carried out on the simplest model of the lattice, a single 
ring of six carbon and six hydrogen atoms, obtaining the equilibrium 
oxygen to lattice distance. These values were then used for the 
calculations on the two ring models directly as well as reoptimized. 
Since the distances, charge distributions, and binding energies did 
net change significantly between the optimized and non— optimized 
calculations, the oxygen to lattice distance was not redetermined 
in the larger structures of the two dimensional lattice models. For 
the three dimensional lattice models, the oxygen-lattice distance 
was again allowed to vary. All calculations were done within a rigid 
lattice approximation. Binding energies for the oxygen atom were 
determined from the difference between the calculated heats of 
formation of the cxygen-graphi te models, the analogous pure graphite 
models and the experimental value of 59.559 kcal/mole for the 
oxygen heat of atomization used in the parameterization procedure 
of MINDO. Electronic charges on the atoms were calculated using the 
MulliJcen population approach. 


54 



Ill « Results and Discussion 

Before the oxygen-graphite calculations were done, the degree 
to which the various models of the lattice describe a perfect 
infinite system was checked. Two quantities often used are the 
valence charge residing on each atom, with a value of 4.00 for the 
ideal case, and the cohesive energy of the lattice, the negative 
of the binding energy with an experimental value of approximately 
5 97/atom. H This latter quantity is of particular interest since 
the previous theoretical studies of oxygen chemisorption on graphite 
resulted in an overestimation of the lattice energy leading to the 
prescription of arbitrarily dividing the binding energies by a 
factor of 5. 

The various clusters chosen to model a single layer of the 
graphite lattice together with the calculated average binding energy 
per atom and the valence charge are presented in figure 5. For the 
three dimensional graphite models with an inter-layer spacing of 
3.35A, the valence charge did not differ from the analogous two 
models. The good agreement between the calculated valence charges 
for the lattice models and the ideal graphite value indicates that 
with proper treatment of edge atoms even smal] clusters can 
adequately simulate the extended lattice. 

Por the oxygen graphite models, results consisting of the atomic 
charges and oxygen binding energies at indicated oxygen to lattice 
distances are summarized in Tables I, n, an d m . The most stable 

position for the oxygen atom is predicted to be above a carbon-carbon 
bond at a distance of 1 «25A for both the independent layer and the 
three dimensional lattice models. The oxygen atom is negative by 
a] most half b C an electron with the nearest carbons in the plane 
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becoming positive. Comparison of the results for the various models 
of thin oxygen-graphite site In Table II], the oxygon-lattice 
distances and charge distributions do not change significantly 
from the poorest racdel (VII] ) where the oxygen Is above nn edge 
bond to the best model (X) where there are at least two shells of 
outer carbons about the bonding site. The local nature of this 
oxygen interaction irlth the nearest two carbons is thus evident. 

The carbon-oxygen distance and charge distributions are typical 
of those found in epoxide molecules. 

Considering the next favored position consisting of an oxygen 
above a carbon atom, in the three dimensional lattice model this 
position gives rise to two sites as described previously. For site a 
where the oxygen is above and below a carbon atom comparj son of the 
calculations with the oxygen at the single plane optimized oxygen 
to lattice distance and at the point midway between the two lnttioe 
layers indicates that the charge on the oxygen (-0.60 electron) 
does not differ significantly from the single plane results. Hcwever, 
consideration of the overall charge distribution summarized in 
Table I major differences between the two and three dimensional 
lattice models arise. 

With the single layer models, the oxygen charge is at the 
expense of the carbon directly below the oxyger. as well as those 
atoms adjacent to this carbon. The presence of the second carbon 
layer in the three dimensional model aliows for the chnrge on the 
oxygen to come from carbons on both layers leading to change c only 
at centers directly bonded to the oxygen and the stabilization of 
the larger oxygen to lattice distance of 1 .67A.. For the single 
layer case, the smaller oxygen lattice distance of 1 .44A is favored. 
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Analogous carbon oxygen Interactions found In molecules are the 
ether compounds with carbon oxygen distances of 1.4 - 1.5 A and 
carbonyls, carbon oxygen distances of 1 ,2A. 

For, the alternate site, site B, where the oxygon Is located above 
one carbon and below the center of the hexagonal ring of carbons 
results indicate a less favored position over that of site A. The 
destabilizing effect of the oxygen unable to interact with the • 
ring atom charge distribution to satisfy bonding requirements is 
evidenced by the positive oxygen to lattice binding energy for the 
three dimensional model (VII) and the two dimensional model (VI). 

Calculations done on models of the remaining oxygen site studied, 
site D, again result in a negative oxygen at the expense of the 
nearest carbon atoms. For the three dimensional model, unlike 
the site A case, placement of the oxygen at the midpoint between 
the lattice layers- does not. result in a lowering of the energy but 
rather an increase indicating a preferred interaction with only a 
single graphite layer. This is further supported by the drastic 
lowering of energy when a bond site position is now made accessible 
to the oxygen from the upper plane as in the cluster XIV. 

Results of the present single lattice calculations compare 
favorably with the three previous theoretical treatments of the 
oxygen graphite system even though different clusters were treated 
as well as the more approximate CNDO 3emi -empirical molecular 
orbital method war; used. The main difference in the three previous 
e;il culatl onu wort: it: the approaches used to mode] the extended 
lattice. Bennett, McCarroll and Messmer? simulated the perfect, 
layer by imposing periodic boundary connections tc United clusters 
of carbon atoms. Dovesi, et aU H used a tight-binding crystalline 
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orbital approach to study periodic overlayer:: of atomic oxygen 
on graphite. The work of Haymj y most closely paralleled the cluster 
approach for the two dimensional lattice models presented .here in 
that hydrogen atoms were used to saturate the dangling bonds at the 
edges of similar clusters. All three calculations were in agreement - 
with the present results shewing the bond site to be most stable 
fell owed by the atom site. Direct comparison of energies was not 
possible since it is well known that the CNDO approximation 
consistently overestimates the binding energy. Additionally, in the 
MINDO approximation parameters had been adjusted to give the best fit 
to experimental enthalpies at 298 *K and hence would somehow incorporate 
macroscopic dynamical effects into a static molecular calculation. 

The closer agreement to the experimental value of the lattice 
cohesive energy using the MINDO approximation compared with the 
CNDO approach may be due to this parameterization. The choice 
of parameters could also explain the discrepency between the present 
results and the calculations of Hayns corresponding to a system 
at C *K predicting oxygen binding at the site D position in contrast 
to the present results. Comparison of the present calculations with 
the previous theoretical treatments is summarized in Tables V and VI. 

Thus from the calculations done on cluster models of oxygen 
ii tercalated graphite using the semi -empirical MINDO/3 molecular 
orbital approaoh predictions of oxygen binding looatlons, oharge 
transfers and relative strengths of bindings have been made. Future 
work on this system would need to consider possible interaction 
between oxygen atoms corresponding to a leas dilute intercalation 
case as well as to allow changes in lattice geometry. An encouraging 
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aspect of the .present study is the apparent ability to model the 
extended lattice using finite clusters once care has been taken 
to adequately saturate the bonding on the edge atoms. Though the 
present results are at best semi-quantitative in nature, they do 
offer information not yet available from experimental data. 
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Table I . Calculated atomic charges for cluster models of the 
oxygen-graphite system, site A. Carbon numbering sequence 
refers to Figure 1. For the three dimensional lattice model, 
the oxygen-lattice distance, D, refers to the lower lattice 
plane. Oxygen binding oxlergy, B, is given in eV. 


I II III IV 


0 

-0.5951 

-0.5971 

-0.5921 

-0.6222 

-0.6042 

-0.5988 

°1 

0.2402 

0.2314 

0.2189 

0.2127 

0.2250 

O .1734 

°2 

0.2273 

0.1644 

0.1617 

0.2336 

0.1358 

0.0820 

°3 

-0.1463 

-0.1325 

-0.1256 

- 0.1018 

-0.1115 

-0.0574 

C 4 

0.1874 

0; 1 669 

0.1620 

0.0978 

0.1281 

0.0655 

°5 

-0.1463 

-0.1560 

-0.1495 

-0.0481 

-0.1288 

-0.0731 

c 6 

0.2273 

0.2602 

0.2556 

0.1422 

0.21 1 1 

0.1351 

C 7 





0.1 1 48 

0.1455 

°8 





0.0566 

0.1262 

°9. 





- 0.0216 

-0.0647 

C 1 0 





0.0296 

O.O 881 

J>U) 

1 .41 39 a 

1 1.4139 

1 .4436 a 

1 .4139 

✓ 

1 .4436 

1 .6750 

E(eV) 

-0.587 

-0.139 

-0.155 

-0.204 

-O .163 

-0.402 

a 0xygen- 

•lattice 

distance optimized. 
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Table II . Calculated atomic charges for cluster models of the 

oxygen-graphite system, site B. Carbon numbering sequence refers 
to Figure 2. For the three dimensional lattice model, the oxygen 
lattice distance, D, refers to the lower lattice plane. Oxygen 
binding energy, E, is reported in eV. All oxygen-lattice distances 
were optimized. 



V 

VI 

VII 

0 

-0.5921 

-0.6556 

-0.5964 

°1 

0.2189 


0.2254 

°2 

0.1617 


0.1543 

°3 

-0.1256 


-0.1252 

°4 

0.1620 


0.1538 

°5 • 

-0.1495 


-0 . 1 474 

°6 

0.2556 


0.2456 

°7 


0.1706 

0.0351 

°8 


0.1706 

0.0357 

°9 


0.1706 

0.0352 

°10 


0.1706 

0.0338 

d(A) 

1 .4436 

1 .0515 

1 .4212 

B(eV) 

-0.155 

2.132 • 

0.565 
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Table III. Calculated atomic charges for cluster models of the 

oxygen-graphite system, site C. Carbon numbering sequence refers 
to figure 3. For the three dimensional lattice model, the 
oxygen lattice distance, D, refers to the lower lattice plane. 
Oxygen binding energy, E, is given in eV. 



VIII 

IX 


X 

XI 


0 

-0.4453 

-0.4577 

-0 . 46 1 8 

-0.4646 

-0.4600 

-0.4606 


0.2817 

0.3063 

0.3033 

0.2838 

0.2993 

0.2989 

°2 

-0.0450 

-0,0638 

-0.0620 

-0.0l6l 

-0.0653 

-0.0651 

°3 

0.0190 

0.0206 

0.0212 

-0.0025 

0.0207 

0.0208 

C 4 

C.0190 

0.0206 

0.0212 

-0.0025 

0.0144 

0.0144 

°5 ' 

-0.0450 

-0.0638 

-0.0638 

-0.0161 

-0.0571 

-0.0568 

°6 

0.2817 

0.3063 

0.3033 

0.2838 

0.2875 

0.2869 

C 7 




-0.0007 



C 8 




-0.0038 



°9 





0.0537 

0.0039 

°10 





0.0326 

0.0031 






0.0047 

0.0051 

°12 





0.0090 

0.0038 

D(A) 

1 .231 8 a 

1 .2318 

1 .251 1 a 

1 .2318 

1 .251 1 

1 .2537' 

E(eV) 

-2.75 6 

-1 .774 

1 

• 

“vj 

CD 

V>l 

-1 .408 

-1 .336 

-1.336 


a Oxygen-lattice distances optimized. 
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Tnble IV . Calculated atomic charges for cluster models of the 

ox yg. Q n-gra phltu system, site D. Carbon numbering sequence refers 
to Figure 4. For the thruo dimensional lattice model, the oxygen 
lattice distance, D, refers to the lower lattice plane. Oxygen 
binding energy, E, is reported in eV. 



XII 

XIII 

XIV 

0 

-0.5302 

-0.5292 

-0.5454 

-0.5198 


0.2001 

0.1972 

0.1248 

0.1270 

°2 

0.3143 

0.2913 

0.1337 

C.1 1 1 1 

°3 

-0.2107 

-0.1958 

-0.0455 

-0.0316 

C 4 

0.2257 

0.1960 

0.0694 

0.0490 

°5 

-0.2107 

0.2913 

-0.0455 

-0.0327 

°6 

0.3143 

0.1972 

0.1337 

0.1097 

°7 


0.0694 

0.1248 

0.1932 

°8 


0.0412 

0.1337 

0.1375 

°9 


-0.0005 

-0.0455 

0.0152 

°10 


0.0128 

0.0694 

-0.0035 



-0.0005 

-0.0455 

0.0331 

°12 


0.0412 

0.1337 

-0.0309 

DU) 

1 .2702 a 

1 .2702 

1 .6750 

1 .6750 

E(eV) 

-0.336 

0.030 

0.362 . 

-0.406 


0 , 

Oxygen-lattice distance optimized. 
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Table V. Comparison of the oxygen atom charge for the oxygen 

atom above a carbon atom, a carbon-carbon bend, and the center 
of the hexagonal ring. All results are for single layer models. 


Present results: 

Atom 

Bond 

Center 

single ring 

-0.60 

-0.45 

-0.66 

two rings 

-0.60 

-0.46 


four rings 

- 0.62 

-0.47 


Reference 7 

-0.40 

-0.25 

-0.19 

Reference 8 

-0.30 

-0.25 

-0.30 
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Table VI . Comparison of the oxygen-graphite equilibrium distance 
for the oxygen atom above a carbon atom, a carbon-carbon 
bond, and the center of the hexagonal ring. All results are 
for single layer models. 


Present results: 

Atom 

Bond 

Center 

single ring 

1 .41 

1 .25 

1 .05 

two rings 

1.44 

1 .25 

' 

Reference 7 

1.50 

1.25 

1 .00 

Reference 8 

1.45 

1 .20 

0.80 

Reference $ 

____ 

1.25 



67 



F 3 guro Captions : 


Figure 1 . Clusters chosen to model oxygen directly above carbon 

atom, site A. For the three dimensional lattice model 'IV 
the upper graphite layer (C7, Cg, Co, C^q ' is Indicated 
by dashed lines. The oxygen atom is above Ci in all 
structures. Hydrogens along the edges have oeen omitted 
for clarity. 

Figure 2 . Clusters modeling oxygen-graphite site B. The oxygen 
atom is above Ci in structures V and VII. The upper 
graphite layer (C7, C8, O9, C10) in structure VII is 
indicated by dashed lines'. 

Figure 3. Clusters modeling the bond site position site C. Upper 
graphite layer (Cg, Cic» Ci 1 , C12) in structure XI is 
indicated by dashed lines. 

Figure 4 . Clusters modeling oxygen position site D. The upper 

graphite layer (C7, C8, Cg, Ciq» dip C^) in structure 
XIII and XIV is indicated by dashed lines with the 
oxygen atom located between the two layers. 

Figure 5 » Clusters modeling the graphite lattice with the calculated 
valence charge indicated. Binding energies are reported 
for each structure in units of eV/atom enclosed in 
parenthesis. 
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1 



XI! 


4 



3 



XIV 


Figure 4 
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3.99 



(- 4 . 72 ) (- 4 , 96 ) (- 5 . 23 ) 

Figure 5 
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